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Abstract 

We study the behavior of a general gravitational action, including quadratic 
terms in the curvature, supplemented by a compact scalar field in 4 + 1 di- 
mensions. The generalized Einstein equation for this system admits solutions 
which are compact in one direction and Poincare invariant in the remaining 
directions. These solutions do not require any fine-tuning of the parameters in 
the action — including the cosmological constant — only that they should satisfy 
some mild inequalities. Some of these inequalities can be expressed in a uni- 
versal form that does not depend on the number of extra compact dimensions 
when the scenario is generalized beyond 4+1 dimensions. 
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1. Introduction. 

The old idea tliat tlie universe might contain more than the observed four space-time 
dimensions has re-emerged recently in novel attempts to explain the weakness of gravity 
compared to the other forces [|l|] and the hierarchy problem , but it was realized earlier 
that such theories might be able to address the cosmological constant problem |^ . The hope 
is that with extra dimensions, the metric might be able, through a non-trivial dependence 
on the extra coordinates, both to accommodate an arbitrary value for the cosmological 
constant and to maintain Poincare invariance in 3 -f 1 of the directions. In their original 
treatment, Rubakov and Shaposhnikov Q found that this idea could be realized with 
two extra dimensions but that the bulk metric contained a singularity. Moreover, their 
system was shown to be unstable by 0. The more recent scenarios involving 3-branes 
reintroduce the cosmological constant problem in the form of a fine-tuning between the 
bulk cosmological constant and the brane tension. The attempts by [^ and to avoid 
this fine-tuning again resulted in naked singularities in the bulk 0. 

We pursue the original idea of Rubakov and Shaposhnikov — that a non-factorizable 
metric that depends on an extra dimension could account for the presence of a cosmological 
constant. In particular, we seek metrics that are periodic in this extra dimension so that 
it can be made naturally compact when the period of the metric is identified with the 
compactification radius. 

The standard Einstein equation for a theory without any 3-branes but with a cosmolog- 
ical constant does not permit metrics that simultaneously satisfy both requirements — that 
the metric is periodic and free of singularities. We therefore consider the effect of including 
higher order terms in the gravitational action which contain two powers of the curvature 
tensor 0, such as a Gauss-Bonnet term |]T^ [|11| |T^. Such terms can be regarded as the 



next natural terms in an effective theory of the gravitational action such as might arise in 
the low-energy expansion of some quantum theory of gravity. 

We show that a general five-dimensional gravitational action including all terms up to 
fourth order in derivatives when supplemented by a compact scalar field permits metrics 
that are periodic in the fifth coordinate and preserve Poincare invariance in the other four 
dimensions. This result does not require any fine-tuning at the level of the parameters 
in the action — including the cosmological constant — other than that they should satisfy 
some mild bounds. Moreover, since these solutions are smooth and contain no singularities, 
either in the metric or in the scalar field, the scenario does not require any 3-branes with 
an accompanying requisite fine-tuning between the brane tension and bulk cosmological 
constant. The only requirement is that the compactification radius should be sufficiently 
small so as not to produce any discrepancies with current experiments. 

In the next section, we review the solution to Einstein's equations in a 5d theory with 
a cosmological constant and a free scalar field. Section three derives the field equations 
from a general action that includes quadratic terms in the curvature. In section four, 
we present several exact solutions of these field equations and further show that with the 
Gauss-Bonnet term alone the metric does not have periodic, smooth, non-singular solutions 
with 4d Poincare invariance. The fifth section shows the requirements that a general R^ 
action should satisfy to admit periodic, smooth, non-singular metrics and discusses a few 
representative cases found numerically. Section six examines a few general properties of 
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including higher order terms in the effective action and discusses theories with more than 
one extra dimension. Section seven concludes. 



2. Background. 

In order to show the importance of higher derivative terms in the action, we first 
examine the solutions for a five- dimensional theory with only the standard Einstein-Hilbert 
action and a free scalar field ,lJ 

S = MI I d^xdy v^ (2A + i?) - ^ j d^xdy ^/^ (Va(/)V"(^) . (2.1) 

Here A and M5 are respectively the cosmological constant and the 5 dimensional Planck 
constant. Qab is the metric for the space-time. We denote the coordinates that correspond 
to the usual space-time dimensions by x^ , where ^, z^, ■ ■ ■ = 0, 1, 2, 3, and the fifth coordi- 
nate by y, with a, 6, c, . . . = 0, 1, 2, 3, y. Ignoring the scalar field, when A > the universe is 
an anti-de Sitter (AdS) space-time while A < corresponds to a de Sitter (dS) space-time. 
We shall often work in units in which M5 = 1. 
When the metric has the form 



ds^ = gab dx'^dx^ = e^^y^rj^^ dx'^dx'^ + dy'^ 
the im and yy components of the Einstein equation for (2.1) are 



f2.2) 



l{A'f-lA" = -K+\{(t)'f 



-i(Ar== -A -!(</>') 

while the equation for the scalar field is 

^" + 2A'^' = 0. 
A solution to equations (2.3)-(2.4) is given by 



/\2 



(2.3) 



(2.4) 



Aiy) 



:^° cos(2^-|A(y-2/o 
(^'(y) = 2v/^sec('2y'^(y 



1/2 



2/0. 



(2.5) 



where Aq and yo are constants of integration, in addition to the trivial </>'(?/) = solution. 
When A < (dS), we are able to obtain periodic solutions. However, the scalar field 
periodically becomes singular and the metric becomes imaginary. In order to remove this 



Our convention for the signature of the metric is ( — ,+,+,+, 
tensor is defined by -R\,^ = daVt, - d^Vl^ + T^.r^ - T^.r^,. 



while the Riemann curvature 



unacceptable behavior from the space-time, we require 3-branes to cut-off the manifold 
in extra dimension before the ill-behaved region is encountered. When A > (AdS), 
only the (t)'{y) = solution is real and the theory can only be compact if we return to 
the original Randall-Sundrum scenario 0. Thus for either sign of A, the scenario must 
contain some 3-branes which necessitates a fine-tuning of the brane tension with the value 
of the cosmological constant. 

In the following, we shall see how the addition of a general R^ action can lead to an 
acceptable dependence of the warp function, A{y), on an extra compact dimension. Note 



that such an action encounters a difficulty when the extra dimension is not compact |13 



(4)/^A^ 



as noted in [|IJ]. If we consider variations about a flat Ad metric, rj^^, — > gliJ{x ), we can 



relate the five-dimensional curvature for (2.2) to the four-dimensional curvature through 

i?(x'^, y) = e-^^y^R^'^Hxf') + ■ ■ ■ , (2.6) 

where R^"^' is the curvature associated with gl^J . By integrating out the extra dimension, 

j d'^xdy^MlR{x^',y)= ( d'^x^~^MlR^^\x^) f dy e^^^^ + ■ ■ ■ , (2.7) 

we determine the effective Planck's constant M4 measured by a four- dimensional observer 
in terms of the 5d Planck's constant M5 through M| = M^ J dy e^^^' . As long as this 
integral is finite we can define a four dimensional effective theory of gravity. This argument 
fails at the next order since 

I d^xdy^R^^ f d^x^/^ (r^^A^ fdyl + ---, (2.8) 

so that in an effective theory, this term receives a correction proportional to the volume of 
the extra dimension. The problem only worsens at higher orders where the [i?^^)]'^ terms 
receive an enhancement of j dy e'^'^~'^^^^y\ One way |T^ to evade this difficulty occurs 
when the term induced in the Ad effective theory is purely topological, as is the case for 
the Ad Gauss-Bonnet term. In this paper we instead consider a scenario with compact 
extra dimensions along with a general R^ action; a Gauss-Bonnet term alone is insufficient 
for our picture, as we show below. 



3. Gravity from a Generalized Action. 

A generic action with up to four derivatives of metric can be written as 

S = S^ + Ml f d^xdy v/^ {2A + R + a R'' + b RabR"^ + c RabcdR'"^"^ + ■■■)• (3-1) 
The additional terms can be interpreted as a squared Weyl tensor, 

abcd'^ — 77 -K — -J^abJ^ -r- Kabcd-tt , [O-^) 
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a Gauss-Bonnet term, 

IT' d2 a p TDdb 1 p ryabcd (o oN 

rj — K — ^JXahJri + Kabcd-ti , l-J-^J 

and a third possible independent term, 

T d2 op -nab , p -nabcd { A\ 

The Weyl term vanishes when the metric is conformaUy flat, as is the case for (2.2), so 
we are free to add some multiple of the Weyl term to (3.4) without affecting the fleld 
equations. Later we shall use /j, and A to denote the coefficients of the third term (3.4) and 
Gauss-Bonnet term respectively: 

fi = 16a + 5b + 4c 3X = 10a + 2b + c. (3.5) 

We include the effect of a scalar fleld, initially through the action for a free fleld, 

S^= f d^xdy v^ (-iA;oV„(^V"(/>) , (3.6) 

although later we shall add a term (Va(/'V"(^)^ which is of the same order as the R^ terms. 
Note that we have included a coefficient ko in equation (3.6). 



Upon varying the action with respect to the metric we obtain |T^ 

~gab^ + Rab ~ 29abR— 29ab [« -R + b RcdR'^ + RcdefR'^ ^ \ 

+2a RRab - 4c RacR\ + 2c RacdfR^''^ - 2(6 + 2c) R^'^Racdb (3.7) 

+ i(4a + b) QabV'R - (2a + 6 + 2c) V„Vbi? + {b + Ac) V^Rab = T„b 

where Tab represents the energy-momentum tensor for the scalar fleld. 

Tab = \ko [Va4>'^b<P - \gab'^c(t)y'(f\ ■ (3.8) 

For a non-factorizable metric of the form of equation (2.2), the sum of the /x/x and the yy 
components of (3.7) yields an equation in which the free scalar fleld does not appear, 

^ \}.A"" + 'iA'A'" + {A"f + AA"{A'f] 

^' ^ (3.9) 

+l\ [A"{A'f + [A'f] - 3(A')' - f A" = -2A, 

so that the difference of these components can be used to determine its behavior, 

kQ{(j)'f = n [A"" -F 2A'A"' -f A{A"f] + ?,\A"{A'f - 3A". (3.10) 

The scalar fleld equation appears in (2.4), but it is not independent being a conse- 
quence of (3.9) and (3.10). When A{y) is a periodic function of y, then (3.10) implies that 
(t)'{y) should also be periodic. Since this equation is non-linear in A{ii), the integral of 
(t)'{y) over one period is in general flnite and non-zero so cj) must itself be compact. Any 
additional dependence of the action on cj), rather than on its derivatives, must be through 
periodic functions. 



4. Analytical Solutions. 



The full fourth-order set of differential equations for a theory with a free scalar field 
and a set of R^ terms with arbitrary coefficients does not admit a simple analytic solution 
except for special cases. In this section, we study two such examples of exact solutions. 
Although these solutions do not produce warp functions A{y) that are smooth, non-singular 
and periodic, they provide partial boundaries for the region of the {A, A, fi} parameter space 
in which we have found such solutions numerically. One of these surfaces requires both the 
fj, and A terms in the field equations. Among this class of solutions is a metric which falls 
off exponentially as y ^ ±00, as in |T^ but without the need of a 3-brane. A second of 
the surfaces, // = 0, represents a theory with only a Gauss-Bonnet term at the R^ order. 
Note that a third boundary lies along the surface A = where we can trivially satisfy the 
field equations (3.9), (3.10) and (2.4) with constant solutions, A{y) = Aq and (^(j/) = (/)o. 
We have found non-trivial periodic solutions only when A < and fi < 0. 



4.1. An Exact Solution. 

We can discover an interesting set of exact solutions by noting that the linear com- 
bination of field equations that eliminates the scalar field, (3.9), does not depend on A{y) 
except through its derivatives and does not contain the fifth coordinate explicitly. Together 
these properties allow (3.9) to be recast as a second order differential equation. 



H 



1^2^2p 1 /rfp\ 



.dP 



di? + 2^lirj +3=P^+^^ + 4P.^ 



+h [Pz^ + z*] - 3z^ -Ip= -2k, 

Zi Zi 



(4.1) 



through the introduction of 



dA 




, , dz 




and 


P{z) - , 


dy 




^ ' dy 



(4.2) 



A simple set of solutions in the full {A, A, n} parameter space is found by substituting 

P{z) =az'^ + b (4.3) 

into (4.1) which imposes the constraints 



A=2- ' 



a 3a -|- 4 



31 1 



^= oT- 



4a + 2 b a + 2 ^ 2 6 a -f 2 

The resulting warp function A{y) is periodic provided that ah > 0: 



(4.4) 



.My) 



cos Vab{y - yo) 



-l/a 



(4.5) 



Note that this solution for A{y), when substituted into (3.10), implies that the scalar field 
is constant, (j)'{y) = 0. 



In general, this solution is not satisfactory since it contains singularities when a > 
and becomes complex for generic values of a < 0. When a = — l/2n, where n is a positive 
integer, the behavior improves so that e"^*-^-* is everywhere real and non-singular when 
A < 0, although for this case e '•^-' vanishes at regular intervals. We shall find in the next 
section that the set of solutions in (4.5) forms a boundary in the {A, A, /u} parameter space 
beyond which periodic, non-vanishing solutions exist when A < and /i < 0. 

We can also use (4.5) to generate a variant of the usual Randall-Sundrum scenario 
T^ , but without including a 3-brane. For example, when a > and A < 0, we find that 



^A{y) _ gAo gg^j^ /'^-|a(a-F2)A {y - yo)) ■ (4.6) 

Although this solution requires fine-tuning the parameters in the action, as given in (4.4), it 
achieves a metric whose zero mode is centered about y = yo and which falls off exponentially 



as y — *> ±cx3. Unlike previous examples ||T^ which used a thick brane to obtain this 



behavior, (4.6) does not require any scalar field. Since this theory has an infinite extra 
dimension, some method of trapping the Standard Model fields near y = yo is further 
required. 

The endpoints of the set of solutions in (4.6), shown in figure 1, have interesting prop- 
erties. Taking a ^ oo (n = 0, AA = ~|): so that the region in which the first derivative 
of the warp function smoothly changes sign has a vanishing thickness, the solution (4.6) 



approaches the warp function studied by Randall and Sundrum [|T^ 



e 



^(y) ^ g^Og-2V^|y-i/o|_ (•47) 



At the other endpoint, a -^ (A = 0, /uA = ■^), the warp function produces a Gaussian 
metric 

eMy)=e^o^^iy-yo)\ (4.8) 

which for A < decreases more rapidly away from y = yo than the standard Randall- 
Sundrum metric. 

4-2. Insufficiency of a Gauss-Bonnet Term Alone. 

A theory with only a Gauss-Bonnet contribution to the R^ action has the advantage 
of producing a set of differential equations with no more than second derivatives of the 
warp function A{y). Unfortunately, for our purpose restricting to such a theory too greatly 
constrains the form of the solutions. Without the /i-terms, equation (3.9) implies that all 
the extrema of A{y), if they exist, are all local maxima or all local minima depending upon 
the sign of A since at the points where A'(y) = 0, A"{y) = ^A. However, the surface 
/U = is still important as it provides a partial boundary to the region in the full {A, A, /i} 
parameter space in which satisfactory, periodic solutions do exist. 

When the metric has the form (2.2), any warp function A{y) that has an extremum 
will encounter a singularity only a finite distance away from it. Integrating the scalar field 
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equation (2.4), (j)'{y) = V— 4Ae~^*^'^^^)~'^°\D and substituting this result into the hnear 
combination of (3.9) and (3.10) that ehminates A"{y) when fj, = yields a first order 
differential equation. Its solution determines A{y) implicitly: 



My - yo) = j^, arctanh (J^) - ^£= arctan (J^ ) (4.10) 



where 



X = Jxl + |AAe-4(^(y)-^o) and xq = \l - |AA. (4.11 



The signs on the right side of (4.10) cannot be chosen separately; however, the minus sign 
should be used to obtain A'{y) = at some point. In the limit A — i> 0, (4.10) reproduces 
(2.5). 

In the interesting A < region, from the expression (4.10) we observe that no accept- 
able solutions exist for any value of A. For example when A < 0, A{y) becomes singular in 
(4.10) at a finite distance from yo- For A > 0, the warp function also contains a singularity, 
but this time it appears in its second derivative. Solving equation (3.9) for A"{y), we note 
that it becomes singular whenever A'{y) ~ ±A~^/^, which again occurs at a finite value of 
y. Appendix A contains a fuller discussion of the properties of the solution (4.10) and the 
location of its singularities. 



5. An Analysis of the Full R^ Action. 

Before describing our approach for finding numerical solutions and presenting several 
representative examples, we summarize with a sketch of the parameter space in figure 1. 
The unshaded area of the figure, which lies in the A < and fx < region of the {A, A, fx} 
parameter space, shows where smooth, non-singular, periodic warp functions exist. We 
can describe another boundary by solving for A as a function of ^ in (4.4), 

AA = -^M±4v/2M-^; (5.1) 

we immediately observe that the shape of the curve of solutions (4.5) is the same for all 
values of the cosmological constant, A < 0, if the A and fx axes are appropriately rescaled. 
This feature corresponds to a rescaling invariance of the field equation (3.9) under y — > ay, 



With (f)'{y) = 0, the only possible solutions are those of the form, 

1/2 



A(y) = ± 



1± a/1- fAA 



y-^ + Ao (4.9) 

vA 



where any of the four possible sign choices is allowed. The case of a scalar field in a A = theory 
was solved in [ p^ . 



/U -^ u'^^, A -^ cr^A and A ^ a ^A, where cr is a real constant. Therefore, we can express 
our results in terms of two dimensionless parameters, {AA, ^K}. 

1 ^ 




Figure 1. A plot of the parameter space {A,/u}. For convenience, we have chosen A = —1 in 
generating this figure. The shaded region does not appear to contain periodic solutions for A{y) 
while in the unshaded region, we have found periodic solutions numerically for arbitrarily chosen 
points. The curve depicts the surface of solutions of equation (4.5); the darker line shows the 
location of the solutions in (4.6). 

A periodic solution for a generic set of values of A, A and // is found by numerically 
integrating the differential equation (3.9). The coordinate y does not explicitly appear in 
the equation (3.9) which moreover only depends on the warp function through its deriva- 
tives. Thus, we can always translate by y ^ y — yo and A(y) -^ ^{v) + ^o to obtain 
another solution. Therefore we can choose A(0) = A'(0) = without any loss of general- 
ity. We also chose A'"(0) = which limits our solutions to those that are even about the 
origin, although we did not find any periodic solutions that are odd when we relaxed this 
constraint. The subsequent evolution of the warp function away from y = then depended 
solely upon the initial value of the second derivative, A"{Q). 

An arbitrary value for A"(0), given some set of values for {AA, /uA}, does not lead to 
a periodic solution. Generically, the warp function tended to reach a singularity at a finite 
distance or to approach asymptotically a solution of the form A'{y) = constant as y -^ oo. 
Between these two extremes, precisely chosen values of A"{0) produce periodic solutions for 
arbitrarily chosen values of A A and /iA within the unshaded region of figure 1. This result 
demonstrates the existence of periodic solutions which naturally select a compactification 
radius without the need to add a 3-brane to the theory or to fine-tune the parameters in 
the action (3.1), as long as they lie within the allowed region of parameter space. 

When the scalar field enters the action only through a free kinetic term, the value of 
ko{(p')'^ is negative for all the periodic solutions that we found numerically. Although this 
result implies that the kinetic energy term has the wrong sign, it appears to be an artifact 
of truncating the scalar action; this unphysical feature disappears when higher order terms 
appear in the scalar action. An example of such a solution is placed at the end of this 
section. 



One region of parameter space that admits a semi-analytic approximation of the solu- 
tion is where A// ~ e^ <^ 1 with A A ^ 1. In this region we can try an oscillating solution 
for the warp function of the form 



A{y) =Ao + e cos{uj{y - yo)) + Aae^ cos(2c^(y - yo)) + 0{e^ 



(5.2) 



for a small amplitude, e. When substituted into (3.9), the only terms that are 0{e) are 
those linear in A{y): 

\lxA"" - |A" = 0(e2) (5.3) 

which, to leading order, requires cu^ ~ ~3//U so that the period of the compact dimension 



IS 



yc ^ 2n^/-ij,/3. 
The 0{e^) terms, which arise from the following terms in (3.9), 



\lxA"" + ?,ixA'A"' + ix{A"f - 3(A 



i\2 



3 /.// 
2^ 



-2A + C(e^ 



determine the value of e in terms of the cosmological constant, 

2 4 A 



(5.4) 



(5.5) 



(5.6) 



and the coefficient of the 0{e^) term in the warp function: A2 = —\. Thus, for the 
linear terms to dominate requires that A^u ~ |e^ <^ 1. Note that since the A terms in the 
field equations contain at least three powers of the warp function, to this order A can be 
arbitrary provided it does not alter the e-expansion: A A ^ 1. 

This semi-analytic estimate is confirmed when we plot the numerical solution to (3.9) 
for A = —1, A = and ^ = —0.01 shown in figure 2. Holding A = — 1 and A = fixed 
and varying n, we have checked that the period and amplitude of the numerical solutions 
is better and better approximated by (5.4) and (5.6) as we let fx ^ 0. 



Aiy) 



0.1- 




Figure 2. A periodic warp function A{y) for A = — 1, A = and /x = —0.01. The initial 
condition is A"{0) = 21.460889. The value of ko{(f)')'^ oscillates about —4, also with a small period 
and amplitude. 

We can find numerical solutions arbitrarily close to the surface /U = 0, for A < | and 
A = —1, and to the surface described by (5.1), provided fi > — ^- Apparently, another 



boundary exists for solutions with fx < — ^ which extends approximately along the line 
/U ~ |A. We have indicated this boundary in figure 1 by extending the shaded region 
beyond the curve (5.1). As A and fj, approach this boundary from above, the numerical 
solutions grow more cusped as in the more extreme example depicted in figure 3 where 
the warp function has broad maxima and sharp, deep minima. The amplitude is also 
significantly larger so that e^^^^ changes by several orders of magnitude. Although the 
figure appears to have cusps at the minima, the function A{y) is actually smooth at these 
points and has no discontinuities in the slope. 



Aiy) 




Figure 3. A periodic warp function A{y) for A = X = fi = —1. The initial condition is 
A"(0) = 5367.89. Despite their cusped appearance, the minima are smooth. 

Most solutions within the unshaded region of figure 1 tend to lie between the extremes 
depicted in figures 2 and 3 as in the example provided by the following sketch. More 
examples are presented in Appendix B. 




Figure 4. An example of a periodic warp function that lies between the extremes represented 
by figures 2 and 3. For this sketch we used A = — 1, A = — 3 and // = — 2. The initial condition is 
^"(0) = 42.2832125. 

Finally, we must address whether the kinetic energy term in the scalar action can have 
the correct sign, /cq = 1, when A{y) is periodic. If we include the additional termtl 



-^A;i(V„(/>V^ 



Other terms of the same order are possible, such as -RafeV^^V (, 
term is the simplest to analyze. 



(5.7) 
but the squared kinetic 
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then the field equations (3.7) require 



/\4 



f, [^A"" + 2A'A"' + ^{A"f + 2A"{A'f] 

+|A [2A"{A'r + {A'r] - u^r - u" = -a + \u<p'f + \k,{<t^' 

^i[A' A"' - \{A"f +2A"{A'f] 

+i\{Ar - liA'f = -A - \uct^'f - ik,{<t^r- 

Since these expressions are only quadratic in ((^'(j/))^, we can solve for 



(5.^ 



{<t^\y))' 



kp 
3/ci 



-1 ± U - 24 






A + IJ.[ A' A'" - -{Ay + 2A!\A!f 



^\\{A!f-\{^f 



1/2' 



(5.9) 



and substitute the result into (5.8) to obtain a differential equation for A{y). This equation 
still admits periodic solutions such as the example sketched in figure 5 for A = —1, A = 0, 
\i — —0.1, /co = 1, and k\ = —0.25 and choosing the minus roota in (5.9). The value of 
ko{4>')'^ is positive so we obtain the standard normalization for the scalar kinetic energy. 



3 + 



2 - 



1 - 






\ A' 



iy) 




Figure 5. A periodic warp function A{y) (solid line) and (j)'{y) (dashed line) for A 
fi = -0.1, ko = 1, and fci = -0.25. The initial condition is A"{0) = 23.77364592. 



-1, A = 0, 



6. Discussion. 

6.1. Higher Order Terms in the Effective Action. 

The gravitational action that we have considered should be regarded as only the first 
few terms of a possibly infinite effective action arranged in powers of derivatives. Therefore, 
we might worry whether higher order terms will spoil the periodic behavior seen in the 
last section. Here we briefly motivate why the existence of smooth, periodic, non-singular 
solutions might be a generic feature of this scenario. 



The choice of the plus root gives a periodic solution that again requires fco < 0. 
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The small A region of the generalized parameter space, formed by the coefficients 
of the terms in the effective action, still should admit warp functions that approximate 
sinusoids as in (5.2), 

A{y) = Ao + e cos{u{y - yo)) + O(e'). (6.1) 

For this approximate solution, the amplitude of the oscillations is assumed to be small, 
(9(e), so terms in the equations of motion with fewer powers of A{y) tend to dominate the 
shape of the warp function. For an action with a general set of B!^ terms included, the 
leading behavior in the e <^ 1 limit is 

fc=3 i 

where the iJik,i are some linear combinations of the coefficients of the R^ terms. When 
the warp function is of the form (6.1), these linear terms have solutions as long as the 
coefficients Hk^i are such that the equation 

^ + i^c.^ + f;$:(-l)VM-^^('=-^^ = (6.3) 

fc=3 i 

has a real root. As before, the cosmological constant can be related to the amplitude e, by 
solving for the 0{e^) terms of the field equations of the full Lagrangian. Thus, when the 
cosmological constant is small, A ~ 0{e^), then the full parameter space contains a region 
in which the warp function is, up to 0{e^) corrections, given by a sinusoid (6.1). As with 
the parameter A in the B? action, the other parameters which appear in the full equations 
of motion, which we denote by \k,i-, only appear in terms that are at least cubic in powers 
of the warp function and so only need to satisfy the weaker requirement that they do not 
ruin the e expansion, \k,i^ ^ 1. 

The problem of establishing the existence of more general periodic solutions becomes 
only more difficult at higher orders. However, well-behaved periodic solutions do seem 
generically to exist for subsets of the terms of the general field equations other than the 
linear terms. As an example, one such subset of the terms in (3.9), 

\^iA"" + [A^ + |A] A"{A'f = 0, (6.4) 

also has periodic solutions when - > — |. One subset of the terms in the R^ action is 

A""" + ci A""{A'f - C2A"{A'Y = 0, (6.5) 

which has periodic solutions for some subregion of the space ci,C2 > 0. This behavior 
also suggests that it is plausible that periodic solutions should exist in some region of the 
enlarged parameter space when higher order terms are included in the gravitational action. 
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6.2. The Scenario in d Dimensions. 

Although for clarity we have concentrated on a scenario in which the universe contains 
3 + 1 infinite dimensions and one extra compact dimension, the picture can be generalized 
to higher dimensions with many of the features found above intact. The most trivial 
modification of the metric (2.2) would be of the form 



ds' 



Aiy) 



d-5 



Tj^jj dx^dx'^ + 2^ dzl 



i=l 



+ dy' 



(6.6) 



for a theory in d total space-time dimensions. Although the warp function A{y) must still 
be a periodic function of y, the Zi directions can be trivially compactified with arbitrary 
compactification radii. The generalized Einstein equation (3.7) for this scenario produces 
the following two independent differential equations, 



H 



I A"" + ^A'A'" + 



3(d-l) 



(A 



//\2 



+ 



(d-lf 



A" {A 



/\2 



+ ~X [A"iA'f + ^{A')^]^ 
-j2 [A' A'' 



d-i ~ 

4 



(d-l)(d-2) / ^/■,2 _ d-2 j^i 

\{A"f + ^A"{A'f] 



d-1 



A (A 



//^2 

/N,4 (d-l)(d-2) 



{A 



i\2 



-A+iA;o( 



-A-iA;o( 



Ll\2 



k'\2 



(6.7) 



As before, only two of the three possible linear combinations of R^ terms contribute. 

d-4 



lj^4{d-l)a + db + 4c A 



[d{d-l)a+{d-l)b + 2c] 



{6.i 



As in the d = 5 case, for a conformally fiat metric such as (6.6) the squared Weyl term 
CabcdC"'^^'^ does not contribute to the equations of motion (6.7). A vanishes for d = A since 
the Gauss-Bonnet term then corresponds to the Euler form and is a topological density. 
Incidentally, for the solution to be periodic in this unphysical case which only has 2-1-1 
dimensional Poincare invariance, /x must be moderately fine-tuned to lie within the interval 

An interesting feature of this more general field equation is that the analytic solutions 
of (4.5) are still solutions when d> b. The generalization of (4.4) is then 

^^^ {d-l){d-2f -^_h{d-l){d-2) 



x4. 



ii = 



i&ji 



36 "" 8 2a + rf 

a(d-2)(3a + d-l) 
h 2a + d-l 
d-2 



(6.9) 



(rf-l)2 h{2a + d-l] 



Note that in (3.5) we have defined A = |A|d=5 to agree with the conventional norniaUzation 
for the Gauss-Bonnet term. 
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Here we have introduced a rescaled set of parameters, {A, A,//}, to emphasize that the 
shape of the curve shown in figure 1 (5.1) is in fact universal — it is the same for arbitrary 
A < and in arbitrary dimension d > 4 when fj, and A are appropriately rescaled: 

9 



XA = --^A±A^J2^A--. (6.10) 

Note that while the location of this set of solutions is independent of d and A when plotted 
in the {AA, //Aj-plane, individual points along this curve do not correspond to the same 
solutions since both the period and the exponent in (4.5) depend on d. 



7. Conclusions. 

A theory with an extra compact dimension and an action with a set of R^ terms and a 
compact scalar field contains sufficient freedom to permit a metric that maintains Poincare 
invariance in 3 + 1 of the dimensions without any fine-tuning of the terms in the action. 
Although the resulting field equations are fourth order, the existence of these metrics can 
be shown numerically. The examples that we have found are smooth, periodic and contain 
no singularities or zeros throughout space-time. As a check, we have confirmed that when 
the numerical solutions for the warp function and scalar field are substituted back into 
the action, the integral over the extra dimensions gives zero, so that the effective four 
dimensional cosmological constant vanishes. We have also found several exact solutions 
which, while they do not themselves provide satisfactory metrics, bound the region of 
parameter space in which exist the desired smooth, non-singular metrics that are periodic 
in the extra dimension. Moreover, these exact solutions continue to exist for an arbitrary 
number of extra dimensions and can be expressed in a universal, dimension-independent 
form. 

As in the original proposal by Rubakov and Shaposhnikov to address the cosmo- 
logical constant problem using extra dimensions, the existence of metrics with periodic 
warp functions that are fiat in the other 3-1-1 dimensions leaves unanswered the question 
of whether they are preferred over other possibilities. For example, we have also found a 
class of metrics that, while still periodic in one dimension, correspond to 3 -|- 1 de Sitter 
or anti-de Sitter spaces in the other components. There is also the question of stability, 
and whether these solutions can be reached from generic initial conditions. The answers 
to these problems are central to the search for a realistic cosmology, but the problems are 
dynamical in nature and a fine turning of fundamental parameters is now not obviously 
required. 
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Appendix A. Singularities in the Solutions for a Gauss-Bonnet Action. 

In section 4.2, we mentioned that the warp function (4.10) for a Gauss-Bonnet action 
encounters singularities after only a finite interval in the extra dimension when A < 0. 
Moreover, the character of these singularities depends upon the sign of A. When A < 0, if 
we let A(y) -^ — oo in (4.10), we discover that this divergence occurs at 



ttVS 1 

±{y-yo) = - 



v/^ 



l + ,/l-|AA + Wl- v'l-lAA 



(A.l) 



Note that the sum on the right side is always real, provided A < and A < 0. 

For A > 0, another type of singularity occurs since then it becomes possible to have 
A"{y) diverge when A'{y) = iA'^/^. (3.9) and (3.10) imply that \O^A'{y) ^ ±^/l±x, so 
inserting x = into (4.10), we see that these singularities also occur at a finite value of y, 



±(y — vo) = +^a/ ^TT arctanh , ^ , _, — h\^^ arctan , ^ ^ . (A. 2) 

Since A!' {\j) = |A when y is an extremum, when A > the warp function only 
contains minima. Thus we have the possibility that A{y) — ^ +00 which, from (4.10), 
occurs as y ^ ±00. However, for AA > |, A"{y) diverges at a finite value for y. 



Appendix B. Several More Examples. 

In the following table, we list the properties of a several more examples of periodic 
solutions. 



A 


fi 


A"{0) 


period Aj 


— A ■ 


i 
10 


i 
10 


8.38996422 


1.2596 


0.5174 


1 


1 

30 


12.20618213 


0.9608 


0.4321 





1 
4 


8.96881784 


2.4476 


1.4186 





1 
10 


8.35597222 


1.2264 


0.4950 





1 
1000 


64.61423 


0.1148 


0.0422 


-10 


-1 


3.7145907 


2.1898 


0.6545 


-1 


-1 


5367.88 


3.005 


7.159 


-10 


-10 


34282 


4.419 


8.501 


-256 


-256 


10439 


9.524 


8.633 



Table 1. A brief list of some of values of A and /x that give periodic, smoothi, non-singular 
solutions when A = — 1. 
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In the table 2, we hold A and ^ fixed while varying A. 



A 


A"{0) 


period Aj 


Tiax ^min 


-20 


494520 


0.7210 


8.0777 


-10 


53680 


0.9502 


7.1593 


-4 


87404 


1.9086 


7.4489 


-4 


47.264364 


1.4166 


1.9249 


-2 


14.5276037 


1.262 


0.7918 


1 

2 


5.32842445 


1.1708 


0.3155 


1 
4 


3.5543619 


1.1586 


0.2167 


1 
10 


2.14671387 


1.1517 


0.1348 



Table 2. A list of values of the cosmological constant A that give periodic, smooth, non-singular 
solutions given X = fi - 



j_ 
w 



Note that a periodic solution for a specific choice of {A, A, n} is not necessarily unique; 
two different values for A"{0) can lead to two different periodic solutions. An example of 
this phenomenon is listed in table 2 and illustrated in figure 5. 



My) 




Figure 6. A pair of periodic warp function A{y) for A = —4 and A = ^ = —0.1. 
condition for the solid curve is A"{0) = 87404 while that of the dashed curve is A"{0) - 
Again, both functions are smooth everywhere. 



The initial 
: 47.264364. 
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